A new fast reconnection model in a collisionless regime 
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Based on the first principles (i.e. (i) by balancing the magnetic field advection with the term con- 
taining electron pressure tensor non-gyrotropic components in the generalised Ohm's law; (ii) using 
the conservation of mass; and (iii) assuming that the weak magnetic field region width, where elec- 
tron meandering motion supports electron pressure tensor off-diagonal (non-gyrotropic) components, 
is of the order of electron Larmor radius;) a simple model of magnetic reconnection in a collisionless 
regime is formulated. The model is general, resembling its collisional Sweet-Parker analogue in that 
it is not specific to any initial configuration e.g. Harris type tearing unstable current sheet, X-point 
collapse or otherwise. In addition to its importance from the fundamental point of view, the col- 
lisionless reconnection model offers a much faster reconnection rate (M c i [ea3 = (c/ui pe ) 2 /{rL.eL)) 
than Sweet-Parker's classical one (M sp = S~ 1/2 ). The width of the diffusion region (current sheet) 
in the collisionless regime is found to be <5 c 'i ess = (c/uj pe ) 2 /r^^, which is independent of global 
reconnection scale L and is only prescribed by micro-physics (electron inertial length, c/ui pe , and 
electron Larmor radius, r£, e ). Amongst other issues, the fastness of the reconnection rate alleviates 
e.g. the problem of interpretation of solar flares by means of reconnection, as for the typical solar 
coronal parameters the obtained collisionless reconnection time can be a few minutes, as opposed to 
Sweet-Parker's equivalent value of < a day. The new theoretical reconnection rate is compared to 
the MRX (Magnetic Reconnection Experiment) device experimental data by Yamada et al. Ji 
et al. 0] and a good agreement is obtained. 

PACS numbers: 52.35.Vd; 96.60.lv; 52.65.Rr; 45.50.Dd; 96.60.pf; 96.60.qe 



I. MOTIVATION 

Magnetic reconnection is the process that enables to 
convert magnetic energy into heat and kinetic energy of 
accelerated particles 3, 4, 5]. In many astrophysical and 
laboratory plasma situations plasma beta (the ratio of 
thermal and magnetic pressures) is much smaller than 
unity. Thus, magnetic reconnection has attracted a con- 
siderable attention as the plasma heating and charged 
particle acceleration mechanism. It is believed that mag- 
netic reconnection is responsible for solar 0, 0] and stel- 
lar flares [1, H| . It has been also observed in the Earth 
geomagnetic tail [T(| ■ 

This work has several motivations: 

(i) The longstanding problems with resistive MHD 
(Magnetohydrodynamic) , i.e. collisional description of 
the magnetic reconnection has been its too slow rate if a 
consistent, the- first-principles approach is used (Til [l2 1: 
or a fast rate [l3[ , but a lack of a proper motivation ([4 j , 
p. 79) (one can still get the fast reconnection if the resis- 
tivity is enhanced in the diffusion region, as one expects 
in physical applications). In the collisionless regime there 
is no known simple, analytical reconnection rate (e.g. 
analogous to Sweet-Parker, or Petschek rates). What 
does exist, however, is a bulk of mostly numerical simu- 
lation work [3 [T! [H, CD! 0. It should be mentioned 
that a significant progress (including analytical) has been 
made in study of the details of collisionless reconnection 
such as the diffusion region structure, outflows, reconnec- 
tion rates [14j, [l9|, |20j, |2l| . However these were either spe- 
cific to a narrow class of collisionless reconnection mod- 
els, namely tearing unstable Harris current sheet; or too 



slow reconnection rates were obtained 12211 . Some numer- 
ical simulation models (Til IT5I \l$L IT7I iTsl [T9I I20I l2Tj| can 
produce as fast reconnection rates as M > 0.5 where M 
is the inflow Alfvenic Mach number. The lack of a simple, 
analytical model of collisionless reconnection (that is not 
specific to any initial configuration) has been somewhat 
hampering the progress in understanding. 

(ii) Quite often, particularly in space plasmas, there is 
a good reason to go beyond resistive, single-fluid MHD. A 
simple line of reasoning can be outlined on an example of 
e.g. solar corona. Fixing coronal temperature at 1.0 x 10 6 
K, Coulomb logarithm at 18.1, the Lundquist number 
(using Spitzer resistivity) is 3.7 x 10 11 . Here L = 10 5 m 
was used. One of the arguments for going beyond resis- 
tive MHD is comparing typical width of a Sweet-Parker 
current sheet S sp = S^ X / 2 L = 0.16 m to the ion inertial 
length. Typical scale associated with the Hall term in the 
generalised Ohm's law at which deviation from electron- 
ion coupled dynamics is observed is, c/u> P i = 7.2 m. Here 
particle density of n — 1.0 x 10 15 m~ 3 is used. Hence, 
the fact that c/ (w P iS sp ) = 44 3> 1 warrants going beyond 
single fluid resistive MHD approximation (a similar con- 
clusion is reached by Yamada et al. [lj in their Figure 12 
based on laboratory plasma experiment known as MRX 
(Magnetic Reconnection Experiment)). 

(iii) Previous results on collisionless reconnection both 
in tearing unstable Harris current sheet [TJ, [T^, [l6|, [T^| 
and stressed X-point collapse [TtI EH has shown that 
magnetic field is frozen into electron fluid and the term 
in the generalised Ohm's law that is responsible for break- 
ing the frozen-in condition is electron pressure tensor off- 
diagonal (non-gyrotropic) components. Thus, a need for 
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inclusion of the electron pressure tensor non-gyrotropic 
components in a model of collisionless reconnection has 
become clear. 

(iv) There is a growing amount of work [T3, EH that 
suggests that in the collisionless regime, on the scales less 
than c/ujpi magnetic field is frozen into the electron fluid 
rather than bulk of plasma. One can write in general 
V e = Vi — j /(en). This relation clearly shows that in 
collisional regime (when the number density n is large), 
the difference between electron and ions speeds dimin- 
ishes V e = Vi = V. However, as one enters collisionless 
regime (when the number density n is small) the devia- 
tion between electron and ion speeds starts to show. In 
Tsiklauri and Haruki [l8[ we proposed a possible expla- 
nation why the reconnection is fast when the Hall term is 
included. Inclusion of the latter means that in the recon- 
nection inflow magnetic field is frozen into electron fluid. 
As it was previously shown in Tsiklauri and Haruki (l7| 
(see their Figs. (7) and (11)) speed of electrons, during 
the reconnection peak time, is at least 4-5 times greater 
than that of ions. This means that electrons can bring 
in / take out the magnetic field attached to them into / 
away from the diffusion region much faster than in the 
case of single fluid MHD which does not distinguish be- 
tween electron-ion dynamics. Thus, a need for inclusion 
magnetic field transport by electrons in a model of colli- 
sionless reconnection has become clear. 

(v) A pioneering work of Yamada et al. [l| has demon- 
strated a clear transition from collisional to collisionless 
reconnection regimes by varying the plasma density, and 
established that in the collisionless regime the reconnec- 
tion rate is much greater than the Sweet-Parker rate. De- 
spite a well motivated explanation of their experimental 
results, one was surprised to see only experimental data 
points and numerical simulation results on their Figure 
12. Thus, we set out with a task of formulating a simple 
model of magnetic reconnection in a collisionless regime. 
In what follows we shall be using well-motivated argu- 
ments of Sweet-Parker model, but shall be applying it to 
the collisionless case. 



II. THE MODEL 

It is instructive to re-iterate key points of the Sweet- 
Parker model. The first step towards the derivation of the 
reconnection rate is that the plasma outflow speed from 
the diffusion region is of the order of the Alfven speed, 
Voutfiow =Va = Ba/y/(i a min = Bo/VMoP- This follows 
from the following consideration: taking the fluid to be 
incompressible and assuming a steady state condition, 
one obtains p. 285 



tflow 



Pi - Po, 



(1) 



where pi and p are thermal pressures inside and outside 
the diffusion region. This pressure difference can be set to 



the magnetic pressure -Bq/(2/i ). From which the above 
result Voutfiow — Va readily follows. 

The second step is applying the continuity equation 
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l = v A s, 



(2) 



where L and 5 are the diffusion region length and width 
respectively. 

The third step is using the generalised Ohm's law (e.g. 
1 p. 108) 

E = —v e x B — + (v e ■ V) v e + r)j, 

n e e e \ at J 

where E and B are electric and magnetic fields, v is 
plasma velocity, P is pressure tensor (3x3 matrix), n 
is plasma number density, m is mass and e is electric 
charge. Subscript e stands for an electron. In Eq.(3) 
we balance two terms advection (v e x B) and resistive 
diffusion (rjj): 



V m fiowB = rjj = r/Bo/inoS), 



(4) 



where j = (V x B)//j,q is used, with V « 1/(5. Note 
that in the (collisional) resistive MHD regime all other 
terms in the generalised Ohm's law are insignificant and 
V e = Vi = V. Finally, substituting (1/5) from Eq.(2) and 
after some simple algebra one arrives at 
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lv A = s- 1 ' 2 , 



(5) 



Where S = [IqEVa/ti, the Lundquist number has been 
introduced. Eq.(5) constitutes the classical Sweet-Parker 
reconnection rate. 

In order to obtain collisionless reconnection rate in 
an analogy with the above derivation of the classical 
Sweet-Parker rate, we now balance the advection (v e x B) 
with electron pressure tensor (V ■ P e / (n e e)), because now 
Spitzer resistivity in the collisionless regime is not large 
enough and also V e ^ Vi = V. Previous results of the 
reconnection in the collisionless regime, both in tearing 
unstable Harris current sheet [ij], [TBI, and stressed X- 
point collapse pjKHj], have shown that: (i) the term in 
the generalised Ohm's law that is responsible for break- 
ing the frozen-in condition is electron pressure tensor 
non-gyrotropic components and (ii) the magnetic field 
is frozen into electron fluid. Kuznetsova et al. [l!| for- 
mulate two dimensional model, which shows that in a 
steady state, at the magnetic null (X-point) where the 
electron flow velocities and magnetic field are zero, the 
only terms in the generalized Ohm's law that can support 
out-of-plane electric field (and hence the reconnection) 
are electron pressure tensor non-gyrotropic components 
(see their Eqs.(4)-(6)): 
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NG 



1 fdP x 



dx 



(6) 



Here the two x- and z-coordinates are along (outflows) 
and across (inflows) the current sheet respectively. Fur- 
ther simplification can be made observing that P xy « P zy 
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but d/dz ^> d/dx because width of the current sheet is 
much smaller than its length: 



E 
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1 P, 



ne 8 C ' ie 



(7) 



where 8 c 'i ess is the width of the current sheet in the col- 



lisionlcss regime {d/dz ~ 
with Eq.(4) we now have 



8 C , 
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Thus in an analogy 



Ey = V e ,inflowBo — 



1 P, 



~y 
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(8) 



Note that the latter equation is similar to Eq.(ll) from 
Hesse et al. [3]. In addition to balancing the differ- 
ent terms, note also that the two crucial differences from 
Eq.(4) are that (i) we use electron speed V e A n flow be- 
cause magnetic field advection on scales less than c/oj p i 
is done by electrons and (ii) width of the diffusion region 
(current sheet) is S c n ess which we shall specify below. 

Let us now specify electron pressure tensor non- 
gyrotropic component P zy . Based in earlier works, which 
use second moment of Vlasov equation (i.e. multiplying 
the Vlasov equation by ViVj and integrating over the ve- 
locity space), Kuznetsova et al. give explicit expres- 
sion for the electron pressure tensor components. Inter- 
mediate expression for P zy in the static case, which is 
generic and yet is not specific to tearing unstable Har- 
ris current sheet configuration, is given by Eq.(14) in 
Kuznetsova et al. 1191: 



P - Pzz dVez 
zv ~ 2n x dz ' 



(9) 



where Q x = (eE>o / 'm e )(z / '8 c n ess ) takes into account linear 
variation of magnetic field with distance in the vicinity 
of the magnetic null. The electron pressure tensor non- 
gyrotropic (off-diagonal) components are generally much 
smaller than gyrotropic (diagonal) ones. The deviations 
from gyrotropic pressure are possible due to electron me- 
andering motion in the regions of weak magnetic field 
close to the x-point [24[ ■ Away from the x-point the par- 
ticles are magnetized and pressure is isotropic. Thus, 
it would be reasonable to assume that electron mean- 
dering motion will be effective up to z ~ n,,e- Here 
ri jE = Vth.e/^ce is the electron Larmor radius. v t h, e — 
^/ksT /m e is electron thermal speed and u> ce = eB/m e 
is electron cyclotron frequency. Further away from the 
magnetic null, when z 3> t\l e then pressure becomes 
isotropic. Thus in the above expression for Q x we set 
z « r L:6 , i.e. tt x = (eB/m e )(r L ^ e /S c >i ess ). In Eq.(9) we 
can also assume the gyrotropic pressure component P zz 
is of the order of the magnetic pressure Bq/(2uq), set 
v ez as the electron inflow speed V et infiow, also estimate 

8/dz W 1/Sc'le 



s . Thus for the P zy we have 
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(10) 



Naturally, the next step is to specify width of the dif- 
fusion region in the collisionless regime, which we do by 



using conservation of mass, but again taking into account 
that now magnetic field is advected by the electrons. 
Thus, in lieu of Eq.(2) we have 



(11) 



where V e: A = i?o/V/ i o nTO e = VA^/mi/nie is the electron 
Alfven speed. 

Thus substituting Eq.(10) into Eq.(8), using Eq.(ll) 
and defining the collisionless inflow Alfvenic Mach num- 
ber, M dless as M c n ess = VeAnflow /V e ,A, after simple 
algebra (also noting that ne 2 /m e = io pe eo and c — 
1/y/fiQeo) we obtain 



(12) 



Note that factor 1/2 has been omitted because this rate 
is a crude estimate. 

Eq.(12) can be regarded as the main result of this work, 
as it provides the collisionless reconnection rate indepen- 
dent of an initial configuration. 

Simple analysis shows that collisionless reconnection 
rate M c n ess = 1.3 x 10~ 4 e.g. for solar coronal parame- 
ters (n = 1.0 x 10 15 m~ 3 , T = 1.0 x 10 6 K, L = 10 5 m, 
B = 0.01T (100 Gauss and hence V A = 6.9 x 10 6 m s" 1 ), 
S = 3.7 x 10 11 ) is two orders of magnitude faster than the 
classical Sweet-Parker rate M sp — 1.6 x 10~ 6 . In the con- 
text of solar flares this means that the collisionless model, 
presented here, effectively alleviates longstanding prob- 
lem of the interpretation of solar flares by means of mag- 
netic reconnection. The reconnection rate is also inter- 
preted as the ratio of Alfven time (ta = L/Va ~ 0.0145 
s) and resistive (or reconnection) times. This means in 
the Sweet-Parker model resistive (or reconnection) time 
is 0.0145/S^ 1 / 2 s = 0.1 days. While our model pro- 
vides reconnection time of 0.0145/(1.3 x 10~ 4 ) s = 2 
minutes, which is commensurate of solar flare observa- 
tions. It should be mentioned that Petschek model also 
can provide an appropriately fast reconnection rate, al- 
though its justification from the fundamental point of 
view is not without a debate 

We now shall compare the master equation Eq.(12) to 
the MRX experimental data of Yamada et al. [Ij ■ Their 
figure 12 presents how the empirical reconnection rate 
scales with the parameter c/ (w P i5 sp ), which is essentially 
the ratio of ion inertial length and the width of the Sweet- 
Parker current sheet. Yamada et al. [l[ have clearly 
established that in the collisionless regime, i.e. when 
c/(u) p id sp ) S> 1 the normalised reconnection rate (after 
starting from the Sweet-Parker rate) attains values much 
larger than unity. In order to compare our scaling law, 
Eq.(12) to the data from Ref.[l|, we need to normalise 
the collisionless reconnection rate by the Sweet-Parker 
rate S 1-1 / 2 . Also, in order for the results to be applicable 
in the collisional regime too, we combine collisional and 

■M c ' less . Noting 



collisionless rates as follows: M = M„. 
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that L = SspS 1 / 2 we obtain, 

M/M sp = 1 + f-^f-V ■ (13) 

TO; 7\L, e \UJpidsp J 

Note that despite a factor m e /rm appearing in Eq.(13), 
M/M sp does not depend on ion mass as l/ui pi oc mj. Also 
the same is evident from the fact that M c q ess (according 
to Eq.(12)) only contains physical quantities pertaining 
to an electron, and M sp is independent of ion mass, so 
should be their sum. Our collisionless rcconnection rate 
seems to fit the experimental data [l[ and the two-fluid 
simulation results [251 ] well in the range c/(aj p iS S p) 3> 1. 
It should be noted that each experimental data point was 
obtained for different plasma parameters. Hence each 
data point has its own S associated with it. S changes 
from roughly 200 to 800 as we move in the lowest mea- 
sured values of c/ (ujpiS S p) = 1.5 to c/(u! p iS sp ) — 10 in the 
horizontal axis (M. Yamada, private communication). In 
order to provide theoretical curves rather than sets of 
the theoretical points we argue that in Eq.(13) the de- 
pendence on S is weak (via cx S sp = LS -1 ' 2 ). Therefore, 
solid curve in Fig.l is produced using a fixed value of 
S = 200. Other parameters used in Fig. 1 were T = 5 
eV, L = 0.4 m and B = 0.05 Tesla. Some deviation in 
the region of c/(w P i<5 sp ) « 1 can be explained by per- 
haps, several factors that were not taken into account in 
our simple model. Let us mention the obvious two: (i) 
Strictly speaking our reconnection rate Eqs.(12) or (13) 
was obtained assuming a steady state, while as Yamada 
et al. [lj mention, the pull magnetic reconnection lasts 
for about 40 /is. Although this is much longer than the 
typical Alfven time < /is, yet one can argue that nearly 
steady state was achieved in the experiment. Thus the 
comparison between the theory and experiment should 
be taken with caution, (ii) Our theoretical reconnection 
rate, naturally, does not include contributions from tur- 
bulence. Despite these shortcomings the match between 
the theory and experiment seems good. 

Based on Eqs.(ll) and (12) it is also straightforward 
to derive the width of the diffusion region (current sheet) 
in the collisionless regime: 

Sc'less = M c ,i ess L — ( ) , (14) 

which is independent of global reconnection scale L and is 
only prescribed by micro-physics (electron inertial length, 
c/ujp e , and electron Larmor radius, r L _ e ). 

Ji et al. [2j have provided MRX laboratory experi- 
mental data and 2D Particle-In-Cell simulation results 
of scaling of the electron diffusion region width, 8 e , with 
the electron inertial length, c/ujp e . It was found that 
the experimental data can be fitted with a straight line 
8c/w pe , while PIC simulation data consistently produced 
a less steep scaling of 1.6c/u> pe . Naturally, we tried to 
apply the diffusion region (current sheet) width formula 
from our model (Eq.(14)) to Ji et al. 0. Useful sketch 
to aid visualisation can be found in Birn and Priest tal, 
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FIG. 1: Solid curve shows collisionless reconnection rate nor- 
malised to the Sweet-Parker rate according to Eq.(13) versus 
c/(uipi8 S p)- Note that Hydrogen, Deuterium and Helium plas- 
mas are all represented by the same solid curve, as M/M sp is 
independent of an ion mass. Data points with error bars cor- 
respond to MRX data Yamada et al. Note that in data we 
do not distinguish between the species. Open square symbols 
are the 2D two-fluid simulation results by Breslau and Jardin 
[25j . Dash-dotted line represents the Sweet-Parker rate. 



p. 91, Fig. 3.1. When plasma inflows in the diffusion region 
ions decouple from electrons at ion inertial length scale of 
w c/ojpi, while electrons deflect from the diffusion region 
(and hence forming it) at electron inertial length scale of 
~ c/ojp e . One of the convincing findings of Ji et al. 0] was 
that the electron diffusion region width, 5 e , is indepen- 
dent of which ion species are used in the experiment (H, 
T>2 or He) . Thus one can conclude that they indeed were 
observing electron diffusion region. This corroborates the 
fact that at the electron inertial length scales the mag- 
netic field transport is done my electrons, as well as the 
currents are carried by electrons. Our width of the diffu- 
sion region according to Eq.(14) is shown as solid curve 
in Fig. 2. A reasonably good fit is obtained, supporting 
the idea that Eq.(14) indeed provides a good theoretical 
expression for the electron diffusion region width. 



III. DISCUSSION 

Let us estimate the coronal heating flux which is pro- 
duced by our collisionless reconnection model. One can 
obtain the heating flux (W m~ 2 ) by taking magnetic en- 
ergy density B 2 /(2/io) (J m -3 ), multiplying it by a typ- 
ical scale L = 10 5 m and dividing by the reconnection 
times which we estimated above, t c n ess = ta/ 'M c n ess = 
114 s (2 min) and t sp = r A /M sp = 8831 s (0.1 days) - 
assuming the reconnection time is a good proxy for the 
energy release time in the solar corona. Fixing B = 0.01 
T, in the case of collisionless reconnection, the heating 
flux is obtained 3.5 x 10 4 W m~ 2 , which is more than 
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FIG. 2: Solid curve shows 5 c 'i esa according Eq.(14). Dashed 
and dotted lines represent 8c/uj pe and 1.6c/oj pe respectively. 
Solid symbols with error bars are MRX experimental mea- 
surements, while open squares are PIC simulation results - 
both from Ji et al. [l. Here T = 5 eV and B = 0.05 Tesla. 



enough to meet the coronal heating requirement; while 
as expected Sweet-Parker model produces an equivalent 
of 450 W m~ 2 . Aschwanden [26|], p. 359 quotes solar 
active region typical heating requirement value of 2000 
W m~ 2 . The obtained collisionless reconnection heating 
flux is an order of magnitude larger than that. Natu- 
rally, neither all of the collisionless reconnection heating 
flux will go into dissipation (by definition collisionless 
process is dissipationless), nor all 100% of the magnetic 
energy will be released by the reconnection process. Even 
if l/10th of it will ultimately end up as the heat, then 
the coronal heating requirement can be met. Thus, the 
importance of the present model for the coronal heating 
seems evident. 

Next, let us comment on the observational conse- 
quences of the collisionless reconnection model in solar 
coronal context. The reconnection time (the flare en- 
ergy release time), t c >i ess = T A /M c q ess cx L 2 (because, 
ta = L/Va, and M c >i esa cx l/L). Thus, if we assume that 
magnetic field and plasma density remains the same as 
we go from small to large flares (this assumption implies 
that for simplicity of the argument Va in t a and c/oj pe 
in Eq.(I2) stay the same), the flare time should increase 
(quadratically) with flare size. There is some ambiguity 
how to define actual spatial size of a flare (Dr. E. Kontar 
of University of Glasgow, private communication): At 
high energies one sees loop footpoints only; Hence is the 
geometric size a distance between the footpoints or total 
(unknown) area? In soft X-rays (6-10 keV) only the loop- 
top source is observed. X-Ray telescope (XRT) on board 
of Japanese space mission Hinode often sees the entire 
flaring loop. However, this mission is quite recent and 
no good flare statistics exists as yet. Temmer et al. [13] 
claim that the flare duration increases with increasing the 
flare importance class, hence with the flare area, as re- 



ported by a number of papers (see their Table 7). Based 
on a statistical analysis of almost 50 000 soft X-ray flares 
observed during the period 1976-2000 Veronig et al. [HI 
show that the scatter plots of the flare duration, rise time 
and decay time as function of the flare fluence (J m~ 2 ) 
show strong correlations. Whether the latter is a good 
proxy of the flare geometrical size is not entirely certain. 
Nagashima and Yokoyama [2!| report a statistical study 
of flares observed with the Soft X-Ray Telescope (SXT) 
on board of Yohkoh in the year 2000. They measured 
physical parameters of 77 flares. Their Fig. (3) plots the 
spatial scale L against the flare timescale and shows that 
the spatial scale L tends to be larger with the increas- 
ing timescale. Thus, in summary, our scaling seems to 
be in agreement with the solar flare observations. One 
should be aware of the fact that in the solar flare ob- 
servations we see only post-reconnection dynamics and 
the geometric size of "particle accelerator", the location 
of actual energy deposition site is never resolved, due to 
poor spatial resolution of the instruments. Some source 
of a concern is also the aspect ratio of the diffusion region 
(i.e. the current layer, where particles arc accelerated). 
From Eq.(I4) it is clear for M c q ess = 1.3 x 10~ 4 the as- 
pect ratio is about 7500 : 1 . The question is whether such 
fine scale structure can survive in the turbulent corona. 
However, since current observations do not have enough 
resolution to prove or disprove existence if such elongated 
current layers, the jury is still out. 

On the laboratory plasma side, comparison of our 
model scaling with the MRX data (see Fig. (2)) needs 
a clarification as to why PIC simulation results (open 
squares) do not fit the experimental data while our sim- 
ple analytical model does. At first sight one might expect 
the opposite to be true because PIC simulation includes 
all relevant physics, i.e. all relevant terms in the gen- 
eralized Ohm's law, time dependence, etc. In author's 
opinion the source of the discrepancy are the boundary 
conditions and the mass ratio mismatch. Ji et al. [2| 
(see their Fig. 3) either use boundary conditions com- 
mensurate to MRX (conducting boundary conditions for 
electromagnetic fields and elastic reflection of particles at 
wall) for unrealistic mass ratio of 400, or open boundary 
conditions for a realistic mass ratio of 1836 (for Hydro- 
gen only). It is no surprise that the outcomes of numer- 
ical simulations depend on correct (appropriate) bound- 
ary conditions used. There are two issues here: first, 
when MRX boundary conditions are used, it may well 
be that a better fit could have been achieved, if the cor- 
rect mass ratio were used (note that experimental data 
points lie on the line (8c/w pe ) that is factor of ~ 5 above 
PIC simulation line 1.6c/w pe . Accidentally such is also 
the mass ratio mismatch 1836. /400 = 4.6. Note, how- 
ever, that the theoretically derived width according to 
Eq.(14) scales as oc l/(ui 2 e r L , e ) cx m e x w ce /v th<e cx 
m e x tti~ x y/ ui e cx y/m e ). Hence if our analytical treat- 
ment is correct, Eq.(14) still does not alleviate the prob- 
lem in full, as -\/1836./400 =2.1 cannot account for the 
factor of 5 difference between the MRX data and PIC 
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simulations. But the trend is in the right direction; Sec- 
ondly, when the correct mass ratio is used (in Ref.0], 
Fig. 3) open boundary conditions are not a good repre- 
sentation of the experiment. Hence the results based on 
open boundary conditions should be discarded. Thus, 
we conjecture that when simultaneously (a) the correct 
mass ratio and (b) boundary conditions commensurate 
to MRX are used, than PIC simulation results will pos- 
sibly follow our theoretical result (Eq.(14)). Until such 
simulation is performed, however, the jury in this issue 
is still out. 

Some clarification is necessary concerning the use of 
notation for L. In Eq.(ll) it plays a role of the electron 
diffusion region length; i.e. the region of space where the 
electron dynamics is decoupled from that of ions. Since in 
the solar flare observations the geometric size of " particle 
accelerator" , the location of actual energy deposition site 
is never resolved, due to poor spatial resolution of the in- 
struments, we use L = 10 5 m, a typical length scale in 
the corona. In context of MRX, we use the length of 
the magnetic scale L = 0.4 m. However, on the MHD 
scale (to which our numerical estimates apply) one also 
expects that at some distance downstream from the X- 
point, electron and ion dynamics becomes coupled again. 
This implies that the outflowing electrons will decelerate 
from the electron Alfven speed to ion Alfven seed. Re- 
cent PIC simulations [3(| El suggest that the electron 
diffusion region length can extend for much longer dis- 
tances downstream than previously thought, and hence 
alleviating, in part, the above problem. 

IV. CONCLUSIONS 

A new model of collisionless rcconnection is for- 
mulated that is based on simple conservation 



laws. The obtained collisionless reconnection rate, 
M c n eS s = (c/ijJpe) 2 l{ r L,eL), naturally does not depend 
on Lundquist number, S, (because it is collisionless) 
and is much faster than the Sweet-Parker rate, but yet 
somewhat slower than the Petschek rate (at least for 
solar coronal plasma parameters). In particular, for the 
same set of parameters e.g. for solar coronal plasma 
the rates are M petschek = 0.04, M c , leas = 1.3 x 10~ 4 , 
M sp — 1.6 x 10 -6 . The main implication for solar flares 
is that if this collisionless rate is used, flare time can be 
as short as a few minutes, that is commensurate with 
the observations. Note that the formulated collisionless 
reconnection model is general and is not specific to any 
initial configuration e.g. Harris type tearing unstable 
current sheet, X-point collapse or otherwise. In differ 
to previous results it relates the reconnection rate to 
simple, generic spatial scales such as electron inertial 
length, Larmor radius and global reconnection length. 
Therefore it is easily applicable to different space or 
laboratory plasma situations. 

Acknowledgments 



Author would like to thank Prof. M. Yamada and his 
team at MRX Princeton Plasma Physics Laboratory for 
kindly providing their experimental data set and for use- 
ful replies to queries. Useful discussions with Dr. G. 
Vekstein and other participants of Fifths International 
Cambridge Workshop on Magnetic Reconnection 2008, 
Bad Honnef, Germany are acknowledged. Author ac- 
knowledges useful discussion of solar flare observational 
aspects with Dr. E. Kontar. Author also would like 
to thank an anonymous referee who helped to improve 
the paper significantly. This research was supported by 
the United Kingdom's Science and Technology Facilities 
Council (STFC). 



[1] M. Yamada, Y. Ren, H. Ji, J. Breslau, S. Gerhardt, 
R. Kulsrud, A. Kuritsyn, Phys. Plasmas 13, 052119 
(2006). 

[2] H. Ji, Y. Ren, M. Yamada, S. Dorfman, W. Daughton, 
and S. P. Gerhardt, Geophys. Res. Lett. 35, 13106 
(2008). 

[3] E. Priest and T. Forbes, Magnetic reconnection: MHD 
theory and applications (Cambridge University Press, 
2000). 

[4] D. Biskamp, Magnetic reconnection in Plasmas (Cam- 
bridge University Press, 2005). 

[5] J. Birn and E. R. Priest, Reconnection of magnetic fields: 
magnetohydrodynamics and collisionless theory and ob- 
servations (Cambridge University Press, 2007). 

[6] D. E. Innes, B. Inhester, W. I. Axford, Nature 386, 811 
(1997). 

[7] T. N. Woods, Adv. Sp. Res. 42, 895 (2008). 

[8] J. H. M. M. Schmitt, F. Reale, C. Liefke, U. Wolter, 

B. Fuhrmeister, A. Reiners, and G. Peres, Astron. As- 

trophys. 481, 799 (2008). 



[9] P. A. Cassak, D. J. Mullan, and M. A. Shay, Astrophys. 

J. 676, L69 (2008). 
[10] A. Runov, W. Baumjohann, R. Nakamura, V. A. Sergeev, 

O. Amm, H. Frey, I. Alexeev, A. N. Fazakerley, C. J. 

Owen, E. Lucek, J. Geophys. Res. 113, 7 (2008). 
[11] P. A. Sweet, Nuovo Cim. Suppl. 8, 188 (1958). 
[12] E. N. Parker, Astrophys. J. Suppl. Ser. 8, 177 (1963). 
[13] H. E. Petschek, in In AAS/NASA Symposium on the 

Physics of Solar Flares, edited by W. N. Hess (NASA 

Press, Washington, DC, 1964), pp. 425-437. 
[14] M. Hesse, K. Schindler, J. Birn, and M. Kuznetsova, 

Phys. Plasmas 6, 1781 (1999). 
[15] J. Birn, J. F. Drake, M. A. Shay, B. N. Rogers, R. E. 

Denton, M. Hesse, M. Kuznetsova, Z. W. Ma, A. Bhat- 

tacharjee, A. Otto, J. Geophys. Res. 106, 3715 (2001). 
[16] P. L. Pritchett, J. Geophys. Res. 106, 3783 (2001). 
[17] D. Tsiklauri and T. Haruki, Phys. Plasmas 14, 112905 

(2007). 

[18] D. Tsiklauri and T. Haruki, Physics of collisionless re- 
connection in a stressed X-point collapse, Phys. Plasmas 



7 



(accepted, inpress) (2008). 
[19] M. M. Kuznetsova, M. Hesse, and D. Winske, J. Geo- 

phys. Res. 103, 199 (1998). 
[20] J. Birn and M. Hesse, Phys. Plasmas 14, 082306 (2007). 
[21] A. Klimas, M. Hesse, and S. Zenitani, Phys. Plasmas 15, 

082102 (2008). 

[22] G. E. Vekstein and E. R. Priest, Phys. Plasmas 2, 3169 
(1995). 

[23] L. Golub and J. M. Pasachoff, The Solar Corona (Cam- 
bridge University Press, UK, 1997). 

[24] R. Horiuchi and T. Sato, Phys. Plasmas 4, 277 (1997). 

[25] J. A. Breslau and S. C. Jardin, Phys. Plasmas 10, 1291 
(2003). 



[26] M. J. Aschwanden, Physics of the Solar Corona an Intro- 
duction (Paxis Publishing Ltd, Chichester, UK, 2006). 

[27] M. Temmer, A. Veronig, A. Hanslmeier, W. Otruba and 
M. Messerotti, Astron. Astrophys. 375, 1049 (2001). 

[28] A. Veronig, M. Temmer, A. Hanslmeier, W. Otruba and 
M. Messerotti, Astron. Astrophys. 382, 1070 (2002). 

[29] K. Nagashima and T. Yokoyama, Astrophys. J. 647, 654 
(2006). 

[30] W. Daughton, J. Scudder, and H. Karimabadi, Phys. 

Plasmas 13, 072101 (2006). 
[31] M. A. Shay, J. F. Drake, and M. Swisdak, Phys. Rev. 

Lett. 99, 155002 (2007). 



